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Abstract 

The collapse transition of an isolated polymer has been modelled by many different ap- 
proaches, including lattice models based on self-avoiding walks and self- avoiding trails. In 
two dimensions, previous simulations of kinetic growth trails, which map to a particular 
temperature of interacting self- avoiding trails, showed markedly different behaviour for what 
was argued to be the collapse transition than that which has been verified for models based 
of self-avoiding walks. On the other hand, it has been argued that kinetic growth trails 
represent a special simulation that does not give the correct picture of the standard equilib- 
rium model. In this work we simulate the standard equilibrium interacting self-avoiding trail 
model on the square lattice up to lengths over 2,000,000 steps and show that the results 
of the kinetic growth simulations are, in fact, entirely in accord with standard simulations 
of the temperature dependent model. In this way we verify that the collapse transition of 
interacting self-avoiding walks and trails are indeed in different universality classes in two 
dimensions. 
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1 Introduction 



Over the past 25 years various lattice models of a single self-interacting polymer chain have 
been analysed in both two and three dimensions (and beyond). These include various types of 
self-interacting self-avoiding walk, self-interacting trails |Tj (lattice paths that can intersect at a 
lattice vertex but not along a lattice edge) and self-interacting random (fully self-intersecting) 
walk models. The fundamental physical phase transition that these models are compared to 
is that of the collapse of single polymer in a poor solvent as the temperature is lowered. The 
question that arises when considering the bulk of these studies is how robust is the universality 
class of the collapse transition. The standard theory [SI 13 El of the collapse transition is based 
on the n — > limit of the magnetic tricritical (/>^ — (j)^ 0(n) field theory and related Edwards 
model with two and three body forces [HlEl) which predicts an upper critical dimension of three 
with subtle scaling behaviour in that dimension. As an analysis of these theories is not exact 
in two dimensions, lattice models form the basis of our knowledge. Analyses of both two- and 
three-dimensional self-interacting trails |H1 EI indicate that they are in a different universality 
class to that of self-interacting self-avoiding walks in those respective dimensions. There is 
no clear understanding of why this is the case if true. A complication of this scenario occurs 
because the numerically most extensive work on trails |Hlini) and that which draws the conclusion 
of separate universality classes, uses so-called "kinetic growth" simulations, or "smart kinetic 
trails", to study the collapse point of self- interacting trails. These simulations are of trails 
produced in such a way that one argues they form a distribution of self-interacting trails at 
one particular temperature. One then further argues from the numerical evidence that this 
temperature is precisely the collapse temperature. Now, it was claimed ^01 that the collapse 
transition associated with "smart kinetic trails" is first-order, though no evidence of this could 
be found in two dimensions. The explanation of this relates to the vanishing of renormalised 
three-body interactions in the smart kinetic trails. It was further suggested in ^01 that studying 
the smart kinetic trails was misleading when considering the full equilibrium self-interacting 
trail model. 

To try to see if these arguments hold in two dimensions we have simulated self-avoiding trails 
over a range of temperatures near the collapse point. The collapse temperature is indeed the 
point onto which the smart kinetic trails map, and we compare these with our earlier results 
concerning two-dimensional trail collapse based on the smart kinetic trail simulations. 

2 The Model 

The model of self-interacting trails (ISAT) on the square lattice is defined as follows. Consider 
all different bond-avoiding paths (^at of length that can be formed on the square lattice with 
one end fixed at a particular site (the set ^n)- Associate an energy —e with each doubly- visited 
site. For each configuration ip]\j count the number m(ipN) of doubly-visited sites of the lattice 
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and give that configuration a Boltzmann weight w™, where u) = exp(/3e). The partition function 
of the ISAT model is then given by 



The reduced and normahzed internal energy 

~ N 



UN = ^Ani) (2.2) 



and reduced and normalized specific heat 
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CN = j^ {{m') - {mf) (2.3) 

are defined in the usual way. 

It was argued in |lHll2j that the smart kinetic growth trails simulates at a specific Boltzmann 
weight uj = uJc = 3 and it was later argued sSi that uJc is indeed the collapse value of the 
Boltzmann weight. 

Let us assume that the phase transition occurring at lOc is critical and furthermore that the 
specific heat diverges at the transition (all previous work supports this second proposition). Now 
we define the exponent a from the divergence of the thermodynamic limit specific heat as usual 
as 

c{uj) = lim CAr~|wc-wr" (2.4) 

Let us return to finite scaling and consider the peak in the specific heat near uJc- Let the peak 
be located at uJc,n and define the width of the transition Au by the difference in the omega 
values at the half-heights (that is, the values of uj that give values of Cm as Cn{'^c,n) /'^)- The 
crossover exponent is defined by 

^^-W^ (2-5) 
where a is a constant, and, in turn, gives us a scaling variable {lo — uJc)N'^ on which to scale for 
data collapse. The standard scaling theory jT^ predicts that the crossover exponent (j) is related 
to the specific heat exponent a via the scaling relation 

2 -a = 1/0 (2.6) 

and that the shift of the peak from its thermodynamic limit value is also governed by (p, that is 

b , , 

|Wc -WcTvl ~ (2.7) 

where 6 is a constant. Moreover, the peak value of the finite size specific heat should then behave 

as [in] 

CNiu^cN) ~ diV-* , (2.8) 

where d is a constant. Using all of the above one can obtain local estimates of the exponents a 
and (j) from this form, and these can be independently verified by studying the width and the 
shift of the transition via the specific heat peak. 
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A further independent test is to study the internal energy as it is predicted to scale as 

Un{cOc,n) ~ ^7ooK) - eA^("-i)* , (2.9) 

where C/qo is the thermodynamic limit internal energy and e is a constant. Usually this would be 
more difficult to analyse as Uooi^^c) is unknown but for this model we can argue that Uoo = 0.4 
as follows from a kinetic growth argument. 

Consider an A^-step loop which occupies M lattice sites. The number of contacts is given 
by the number m of doubly visited sites, and we have m = N — M. Any site of this loop could 
have been the starting point, and in order for a site to be visited twice, the loop must not have 
been closed at the first return visit. The probability of not closing upon the first return is 2/3, 
so that for large loops m/M 2/3, from whence it follows that m/N 2/5. 

In our previous work |B] we estimated 

= csst":"^ (2.10) 

and so a ~ 0.86. This is equivalent to the value of a(j) ~ 0.76 for the exponent describing the 
divergence of the specific heat peak. We note that the established values for self-interacting self- 
avoiding walks in two dimensions ^Ij, that is, (f) = 3/7 and a = —1/3, imply that the specific 
heat does not diverge! On the other hand, if there was a first order transition a = (p = 1. 

3 Results 

We have used the now standard PERM algorithm |T5] to simulate ISAT at various fixed tem- 
peratures around = 3 for trail lengths from 2^^ = 32768 (denoted as 32K) in factors of 4 up 
to 2^1 = 2097152 (denoted as 2048K). We chose the range of temperatures by first simulating 
at u; = 3 and then reweighting the obtained histogram to give an estimate of the location of the 
specific heat peak uJc,n- We then repeated the simulation at uJc,n- The multi- histogram method 
|16j was then used to give data throughout the transition region. The simulations at = 3 gave 
almost identical results to those additionally using the peak data. The quantities of interest 
such as Auj and Wc,Af were re-estimated. 

We attempted to scale all the data using a consistent set of exponents. The best fit for the 
specific heat data was obtained for 

a4> = 0.66 so a = 0.80 and (p = 0.83 (3.1) 

while for the internal energy data the best fit used 

(j) = 0.85 so a « 0.82 and acj) = 0.7 (3.2) 

Clearly there is an error of at least 0.025 in these estimates and the statistical spread confirmed 
an error of about 0.03. We have used the scaling relation (|2.6|) for consistency. The first figures 
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show the the scahng of various quantities using those exponent assumptions. In figure (CJ the 
scaled difference of the internal energy to its thermodynamic value (C/qo — U]\f)N^^~'^^'^ and scaled 
specific heat CnN~°"^ are plotted against the scaling variable — uJc)N'^ . These demonstrate 
consistency in three ways: between the specific heat, the internal energy and the horizontal 
scaling variable all at once. 

We have estimated the crossover exponent independently via both the shift {uc^n — ^cj and 
the width Aw of the transition. Figure (jSJ shows plots of the appropriate scaling combinations for 
these two quantities against using the estimate (p = 0.84 (in between our previous estimates), 
being the best value using these plots. 

Finally, in Figure © we demonstrate the convergence of the scaling by plotting scaling 
combinations that should be constant in the absence of corrections-to-scaling. These give us 
confidence in the results from the other figures. 

We therefore conclude that our new simulations away from the "smart kinetic growth point" 
are in agreement with the exponent estimates of the smart kinetic growth simulations with 
perhaps slightly smaller values of (j) and a as previous central estimates. We now estimate 

a = 0.81(3) and = 0.84(3) (3.3) 

4 Conclusion 

We have simulated self-interacting self-avoiding trails on the square lattice up to lengths of 
2, 097, 152 for a range of temperatures around the collapse transition temperature. We conclude 
that the results are in good agreement with earlier simulations based on smart kinetic trails. 
They demonstrate a phase transition which is neither of the type displayed by self-interacting 
self-avoiding walks nor is it first order, although it is very strong, as predicted by an analysis of 
the smart kinetic growth trails. 
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Figure 1: Scaled internal energy gap (U^ — Un)N^^~°'^'^ and specific heat CnN~°"^ versus scaled 
{oj — uJc)N'^ for lengths 32K, 128K, 512K, and 2048K, comparing the individual data points 
with the results from the multi-histogram method. 
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Figure 2: Scaling of the transition: shift and width of the collapse region. Shown are the scaling 
combinations N'^{lJc^n — "^c) 0''nd N'^Alo versus N. We have that Wc = 3 and have used the 
estimate of (p = 0.84. 
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Figure 3: Scaling of the transition: height of the specific heat peak. Shown are the scaling 
combinations CN{i^c,N — ^c)'^ o-nd CnAuP' versus N . We have that lJc = 3. 
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